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Abstract 

We study the Toulouse limit of the multi-channel Kondo model defined 
as the limit of maximal anisotropy which can be achieved without changing 
the infrared behaviour of the model. It is shown that when the number of 
channels 'k' exceeds two, the interactions do not vanish and the Toulouse 
limit is described by a non-trivial field theory. Considerable simplifications 
take place only at k = 2 where the Bethe ansatz reproduces the results 
obtained by Emery and Kivelson [Phys. Rev. B 47, 10812 (1992)]. 

PACS numbers: 72.10.Fk, 72.15.Qm, 73.20. Dx 



1 Introduction 

In his original paper Toulouse suggested that the anisotropic spin- 1/2 Kondo 
model greatly simplifies at some special value of anisotropy, where it becomes a 
model of free fermions 1 . Later it was discovered that a similar point exists for 
the two-channel Kondo model, but the limiting theory is a theory of free Majo- 
rana fermions 2,3 . As we shall demonstrate below, if the number of channels 
exceeds two, the Toulouse limit does not correspond to a theory of free particles. 
For example, it has been established by Fabrizio and Gogolin 4 that for k = 4 
the Toulouse limit describes the boundary sine-Gordon model with the scaling 
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dimension of the boundary term A = 3/4. The latter is a non-trivial field theory 
whose infrared behaviour is governed by the same scaling laws as for the isotropic 
4-channel Kondo model 5 . Under these circumstances one may wonder how to de- 
fine the Toulouse limit. We think that the logical choice is to define it as the 
limit of maximal anisotropy for which the infrared behaviour still remains the 
same as for the isotropic model. Since the Toulouse limit plays a considerable 
role in many studies of the Kondo models (see Refs. 2,3 and 6), we undertake to 
study it in detail using the Bethe ansatz. 



2 Derivation of the Thermodynamic Equations 

The Bethe ansatz equations for an anisotropic k-channel Kondo model are similar 
to those for the Heisenberg chain with spin k/2 and are given by (see, for example, 
Refs. 6-8): 

M 

[e k (r]] u a )] N e 2 s(m u a ~ 1/flO = II ^(Vi U a ~ u b) (1) 

6=1 
M x 

E = Y,^e k ( V ;u a ) (2) 

a=l 

sinhlWu — in)} 
e n { m u)= , ' 3 
smh[r]{u + in) J 

where S is the impurity spin, g is the Kondo coupling constant, rj is the anisotropy 
and N is the length of the system and M = kN/2 + S — S z is the number 
of up spins. The universal relationship beetween the quantities g and r\ and the 
parameters of the Hamiltonian exists only in the limit of weak anisotropy r\ << 1. 
As we have mentioned above, the difficulty in determining the Toulouse limit is 
resolved if we define it as corresponding to the maximal value of r\ at which the 
IR fixed point of the Kondo model still belongs to the same universality class as 
at rj — > 0. The periodicity of the trigonometric functions in Eqs.([l|) suggests that 
for k > 1 the Toulouse limit is realized at r] — n/2(k + l/u) {v — > oo). Here 
we percieve a profound difference between the single-channel (k = 1) and multi- 
channel Kondo models. In the single-channel case the Toulouse limit is defined as 
the point where the scattering vanishes. It occurs at rj = it /A where the right hand 
side of Eqs. (|T|) becomes 1. At k > 2 the point n = 7r/4 cannot be continuously 
connected with the isotropic point rj — > 0, since at rj — n/2k the left hand side 
becomes equal to 1, which marks a discontinuity. Therefore the Toulouse limit 
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defined as the limit of maximal anisotropy is no longer equivalent to free particles 
and corresponds to a non-trivial field theory. Thus, as we have already mentioned, 
according to the arguments given in Ref. 4, at k = 4 the Toulouse limit describes 
the boundary sine-Gordon model with the scaling dimension of the boundary 
term A = 3/4. 

The Toulouse limit has been considered in Refs. 7-9 in the context of the 
model of the Anisotropic Principal Chiral Field. The most detailed derivation 
of the thermodynamic Bethe ansatz (TBA) equations is given in Ref. 9. The 
equations for the Kondo model are essentially the same as for the Principal Chiral 
Field, but, nevertheless, some details need explanation which is presented in this 
article. We will show that the TBA equations in the Toulouse limit have the 
following form: 

e n = s*ln(l + e e "- 1 )(l + e e " +1 ) 
+Sn tk - 2 s * In (1 + e efc ) + <5 nifc _is * f(v), n = l, ...k — 1 (4) 
e k = s * In (1 + e efc ~ 2 ) - 2 exp(-7rv/2) - s * f(v) (5) 

/°° 2 r i 

dvs(v H — \nT K /T) In 1 + e t2s{v) (6) 
-oo 7T L J 



where 



/oo 
du 
-oo 
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4cosh[7r(f — u)/2] ' 

the function f{y) is defined later in the text (see Eq.(|2"9|) and the Appendix) and 
the Kondo temperature T K is defined as: 

7rA 

T K = lim — exp(-7r/2#) (7) 

At zero magnetic field f(v) = 0. When taking the Toulouse limit one has to be 
careful to keep the Kondo energy scale, T K , finite. This is because at v — ► oo the 
left hand side of Eqs.(|T|) formally becomes equal to one. In order to get non-trivial 
solutions one has to renormalize the cut-off, A ~ v. 

As it was shown by Takahashi and Suzuki 10 for k = 1, and later for general k 
by Kirillov and Reshetikhin 8 , at rj = n j{k + Eq.(|l]) has the following string 
solutions: (i) there are k — 1 strings of lengths n — 1, 2, ...k — 1 whose energies we 
denote as e n , (ii) there are kink and antikink excitations E a , o = ±, (iii) there are 
v — 1 breathers with energies Kj, j = 1, 2, ...v— 1. The total number of excitation 
branches is equal to k + v. Kinks and antikinks have spin 1/2 with a Landee 
factor which is a function of the interactions, whilst the breathers are spinless. 
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The TBA equations for the strings are: 

e n = s*ln(l + e £ "- 1 )(l + e e " +lN 



where 



s * 



.-E 



3 

cosh[(l — jj u)uj] 



, n = 1, ...k — 1 



* J ' cosh uj 

In the limit v —>■ oo this kernel becomes a delta-function. 

The TBA equations for the breathers are: 



In (1 + e*') - B jm * In (1 + e" K " 



2v 



— sm(TTj/2u)e- nv/2 + bj*J2 ln i 1 + e ~ £<T ) - & * In (1 + e 6 *- 1 ) 

where the Fourier transforms of the kernels are given by: 

„ , . , , . . cosh[(l — jlv)uj\ smhimui I v) . . 
B jm (u) = 2coth(cu/z/) ^ A , " — » > ™) 



(9) 



cosher 



Bjm B rn j 



bj(uj) 



coth(u; / v) sinh(ju; / v) 
cosh a; 



In the limit i/^oowe have 

B jm (v) = 2min(j,m)5(v), bj(v) = 2js(v) 



(10) 



In this limit Eqs.(^) can be solved and the functions £j expressed in terms of 
other excitation energies. 

The TBA equations for kink and anti-kink excitations are 



: e -™/2 _ 

7T 



g L <?H/T - ]T s * ln (1 + e ej ) + J2 h i * ln f 1 + e ~ K 



;i2) 



where 



K(io) 



sinh[(l — u^^uo] 
2 cosh u sinh(a; / v) 
The free energy has the following form: 



bulk 



-jT 2 I dn-.— ' 2 



7T 



^(l + e- Bff ) +z/~ 1 ]Tjm(i + e _ "' 



d^[^ + -ln(T^/T)]ln (l + e £2s( ^ 

-oo 7T v 



(13) 
(14) 
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where A = N/L and L is the length of the system. The Kondo temperature Tk 
is defined by Eq.(|7|). 

Now we shall consider the limit v —>■ oo. As we have said above, we shall 
keep Tk finite when taking this limit, which implies keeping the combination 
A/u finite. Let us consider the equations for E a first. At v —>■ oo the kernel 
K{uj) — > vK[lS) where 

~ . tanhtu 

K(u>) = 

V ; 2u 

Since the right hand side of Eq. (|l2D is proportional to u, E_(v) goes to infinity and 
only E + (v) remains finite and only in a finite magnetic field. The corresponding 
equation is 

K * In (1 + e - E+{v) ) = ~g L H/2T - - e ~ nv/2 (15) 

where 

ql = lim — (16) 

u^oo y 

The latter quantity is the effective impurity Landee factor which is difficult to 
determine in the Bethe ansatz framework. Similar difficulties are common for 
all models with linear spectrum and originate from the fact that in the Bethe 
ansatz one chooses as the reference state a ferromagnetic state with a very high 
energy (see the discussions in Ref. 11). According to Refs. 3 and 6 the impurity 
magnetic susceptibility vanishes in the Toulouse limit, which indicates that cjl = 
0. However, this result was obtained under the condition that the impurity and 
the conduction electrons have the same Landee factors. For the general case 
9im P 7^ 9c, where g imp and g c are the bare Landee factors of the impurity and the 
conduction electrons, it was demonstrated by Fabrizio et al. 12 that 

9l = gimp - g c (17) 

It make sense to consider gi mp ^ g c and keep c/l finite. In what follows we shall 
denote f(v) = In (l + e -JS +W) . The solution of Eq.flTSD is given in the Appendix. 
The v — > oo limit of the TBA equations for breathers are given by 



In (l + e^) -2£ min(j, m) In (l + e" KmW ^ 

m 

j (V 7 ™/ 2 + 2s * /(«)) - In (l + e ek - l{v) ) {11 
Inverting the kernel we transform these equations to the more familiar form: 



Kj(v) = - In (l + e ^- lM ) (l + e K >+ l( - v) ) - -S j;1 In (l + 



lim iM. = ( e -W2 + 2s * f( v )) = e (v) (19) 
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The general solution of this system was obtained by Takahashi 13 : 

f sinh[±e n (j + 1) - u(v)} ) 2 

1 + e Kj 



I 



sinh(ie ) 



(20) 



In our case the function u(v) is determined by the boundary condition at j — 1: 



1 + e 



sinh(|e ) 
sinh(ie - u) 



(21) 



Using this result we can calculate the infinite sums appearing in other TBA 
equations: 



^ln(l + e K = In 

j 

£im(l+e^ 



-eo/2 



sinh(e — it) 



sinh(|e — u) 
') = -ln(l -e 2u ' t0 

ln(l + e tk )/(l -e t0 



(22) 



(23) 



Here we introduced the new energy e& in order to get rid of u. We shall also 
redefine eu-i introducing a new function d k _ x : 



l + e e *-i = (1 + e efe ) (l + e 6 * 



(24) 



Substituting Eqs.(^),(^) and Eq.(|2§ into Eqs.(9) we obtain Eqs.(§), and (§ 
with /(f) given by Eq.([L^). 



3 Thermodynamic Properties in Zero Mag- 
netic Field 

The system of TBA equations (|5p and (|5p really have different topologies at k = 2 
and k > 2. In the former case the equations decouple and we have 

e 2(v) = —2exp(—nv/2) — s * f{y) 

ei = s* f(v) (25) 

The function / is independent of ei >2 , the explicit expression for it is derived 
in the Appendix (see Eq.(38)). Substituting Eq.(38) into Eq.(6) we obtain the 
following expression for the impurity free energy: 

F imp = -T J dvs[v + -lnT K ) In jl + exp[s * f(v + - InT)]} 
= ~ jH jf^ In {l + expi^(~g L H/T K y + x* - *)]} (26) 
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This formula describes the free energy of a Majorana fermion with dispersion 
E(v). Since / ~ H, E(v) vanishes in zero magnetic field together with the 
impurity part of the specific heat. Thus we reproduce the results obtained in 
Refs. 2,3 and 6. 

At k > 2 no simplifications take place. Neither the impurity contribution to 
the specific heat disappears at H = 0, and the Toulouse limit is described by a 
non-trivial field theory. We can solve the TBA equations only approximately for 
T » T K and T « T K . At T << T K we can expand Eqs. (|]) and (^) in powers 
of gk = In [1 + exp(efc)] and /, as was done in Ref. 14. The answer is: 

$(Jfe-l) 
- 2)1(25) 



s ^ 1+e e 2S (,-^ln(T K /T)) 

x [$(2S - l)/ 25 +2 - $(2S + l)/ 25 ] * 



+ H // 2 ! * f( V ) 

4 cos / (7t/k + 2) 



, . . sinhna; 
L(uj) = . (27) 
smn(fc + 2)uj 

where 

= sin [7r(n + l)/(fc + 2)] 
1 J " sin[7r/(fc + 2)] 

The first non- vanishing pole of the kernel is at u = —2m/{k + 2) which gives 

the correct low temperature asymptotics for the impurity free energy: 

F jmp = -ATiT/Tx) 1 '^ 2 ' 
^/ d „e-~ [9tW + _J^L_„ /W] (28) 

4 Thermodynamic Properties at T = 

/<">(«) + /°° if (« - u)/«(tt)d« = ^ - -e~™ /2 (29) 

JQ 2i 7T 

The solution of this equation is obtained in the Appendix. Substituting Eq . (|38|) 
into the expression for the impurity free energy (|2q) taken at T — > and dif- 
ferentiating the result with respect to the magnetic field, we get the following 
expression for the impurity magnetic moment: 

Mm = f dxsin(2S7r/fc) qlH 

vr J cosily - cos(25tt//c) JglH 2 + T^exp(-7ikx) 
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For k = 2, S — 1/2 this integral can be calculated: 



imp 



h 



In (V 1 + Vh- 2 - l) 0(1 - /i) + sin" 1 y/l - h- 2 6(h - 1)] (31) 



irj\l - h 2 



For k = A, S = 1 the integral for the magnetization also looks rather simple, but, 
nevertheless, cannot be expressed in elementary functions: 



In both cases h = g^H/Tx- 

5 Conclusion 

Here we just give a brief summary of the picture described in the main text. It 
turns out that the Toulouse limit is very different for k = 2 and k > 2. In the 
former case one can solve the TBA equations explicitly for any temperatures and 
magnetic fields. It is also true that at H = the impurity contribution to the 
free energy becomes trivial Fi mp = — T In \/2- This result holds for any H when 
Qimp = 9c (see Refs. 3, 6). Nothing of that sort happens at k > 2. The TBA 
equations remain non-trivial and no explicit analytic solution is available. 
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7 Appendix 

Eq. (|l5D is quite special in that respect that its solution can be expressed in 
elementary functions. Namely, one can treat this equation as a limiting case of 
the equation 




(32) 



Q 



/( ~ 

duln | coth[— (f — u)]\f(u) 



irgiH + Am cosh(7ru/2) 



(33) 



-Q 
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The latter equation becomes Eq.(|i~5|) at m — > 0, mexp(nQ/2) =const with 

f(v) = Tf(v + -\nT) 

7T 

Differentiating this equation with respect to v and using the substitution 
f(u) = [cosh(7r?i/2)]" 1 F(M), x = tanh(7ra/2) 

we reduce it to the canonical form 

rB dyF(y) x 

P / = g(x); g{x) = -nm- - (34) 

J-b x — y 1 — x A 

This equation has the following general solution 15 : 

F(«) = -\JW^P j B , , (35) 

7T 2 J -b yy — x)\J B l — y l 

In our case we get the following expression for f(u): 

\B 2 cosh 2 (7rn/2) - sinh 2 (7ra/2)l 1/2 

/(«) = (1-^)1/2 ( 36 ) 

The limit 5 can be determined from the original equation. In particular, substi- 
tuting there the solution (B6[) and putting v = we obtain the following condition: 



1/2 



/ dw In | coth(vn;/4)| [sinh 2 (7rQ/2) cosh 2 (vn;/2) - cosh 2 (vrQ/2) mnh 2 (7ru/2) 
J-Q 

= -4 + 7tHg L /m{37) 



where 5 = tanh(7rQ/2). In the limit H/m — > oo when Q — > oo we obtain from 
Eq.(0): 

me Q7r/2 = 

and 

s * / = M (_ e -W2 + V ^^J) (38) 
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